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Abstract 

In the ferromagnetic phase of the q-state Potts model, switching on an external magnetic field in¬ 
duces confinement of the domain wall excitations. For the Ising model {q = 2) the spectrum consists of 
kink-antikink states which are the analogues of mesonic states in QCD, while for q = 3, depending on 
the sign of the field, the spectrum may also contain three-kink bound states which are the analogues 
of the baryons. In recent years the resulting “hadron” spectrum was described using several different 
approaches, such as quantum mechanics in the confining linear potential, WKB methods and also the 
Bethe-Salpeter equation. Here we compare the available predictions to numerical results from renor¬ 
malization group improved truncated conformal space approach (RG-TCSA). While mesonic states in 
the Ising model have already been considered in a different truncated Hamiltonian approach, this is 
the first time that a precision numerical study is performed for the 3-state Potts model. We find that 
the semiclassical approach provides a very accurate description for the mesonic spectrum in all the 
parameter regime for weak magnetic field, while the low-energy expansion from the Bethe-Salpeter 
equation is only valid for very weak fields where it gives a slight improvement over the semiclassical 
results. In addition, we confirm the validity of the recent predictions for the baryon spectrum obtained 
from solving the quantum mechanical three-body problem. 


1 Introduction 

In this work we consider the scaling field theory corresponding to the (7-state Potts model in the vicinity 
of the critical point separating the paramagnetic and ferromagnetic phases. In the absence of external 
magnetic field, the field theory in the ferromagnetic phase has a spectrnm consisting of kinks (domain 
wall excitations) and is integrable with a known factorized S matrix [1]. Switching on a weak external 
magnetic field indnces a linear potential between the kinks, leading to their confinement. 

In the case of the Ising model {q = 2) this scenario was proposed in [2], and has been investigated 
nsing several methods since then. One approach consists of a numerical Hamiltonian truncation method 
[3], which allows the numerical determination of the resulting meson spectrum. On the other hand, there 
are several approaches allowing predictions for the spectrum, with the first prediction already made in 
[2], which is in fact equivalent to finding the spectrum of quantum mechanical bound states in a linear 
potential. Form factor perturbation theory which treats the strength of integrability breaking (such as 
introduced by the magnetic field) as a small parameter can also be used to establish confinement [H [5]. 
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The meson spectrum can also be computed using the WKB method |6], and a low-energy expansion can 
be obtained from the Bethe-Salpeter equation la El El El |9]. The validity of these approaches and their 
consistency has been established not only on the qualitative, but also at a quantitative level to high 
precision. 

In the q = 3 case the spectrum is expected to be richer and contain baryonic three-kink bound states 
besides the mesonic kink-antikink ones |inj . For the 3-state Potts model, the Hamiltonian truncation 
method used in mu is not available, since it requires the knowledge of exact finite volume form factors 
of the magnetization operator in the off-critical theory in zero magnetic field. However, an alternative is 
provided by the truncated conformal space approach (TCSA) introduced in |llj . Since construction of 
low-energy spectrum in the TCSA does not depend on the assumption of integrability, nor on the existence 
of small parameters, it can be readily applied to non-integrable models, as performed for the q = 2 and 
3 cases in [12] and m, respectively; for more recent applications to other non-integrable quantum field 
theories cf. El I15| . In the case of the 3-state Potts model, a previous study has already confirmed 
qualitatively the expected mesonic and baryonic spectra |13) . 

More recently precise quantum mechanical and semiclassical predictions for the mass spectrum have 
been obtained |16[ EZl, which makes it possible to investigate the spectra in more detail. For this pur¬ 
pose a more advanced implementation of TCSA is needed, which is based on the renormalization group 
improvements recently introduced and developed in |18[ [T9l [20l 121112^ [23] . In this work we discuss the ap¬ 
plication of this RG-TCSA method to the Ising and 3-state Potts models and compare it to the theoretical 
predictions. 

The outline of the present paper is as follows. In Section]^ we review the known results about confine¬ 
ment in the scaling g-state Potts field theory. In Section we briefly describe our TCSA implementation. 
Section [^contains our results, starting with the Ising case as a testing ground for our numerical procedures 
and then turning to the 3-state Potts model. In Section we draw our conclusions. 


2 Confinement in g'-state Potts field theory 

2.1 Brief overview 


In this section we briefly overview the phenomenology of the g-state model based on [10] • The g-state 
Potts field theory is the scaling field theory of the g-state Potts model which is the generalization of the 
Ising model with q different values (colours) of the lattice variables [23]. The lattice Hamiltonian can be 
written as 
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1 

f 




s{x),s(y) 


H ^ ^s(x],q 


( 2 . 1 ) 


{^,y) 


The first term is the standard nearest-neighbour interaction and second term is a generalized magnetic 
field in the “direction” of the g-th colour. Without magnetic field {H = 0) the theory has Sq permutation 
symmetry and has a critical temperature T = Tc, below which the system is in the ordered (ferromagnetic) 
phase while above the system is in the disordered (paramagnetic) phase. 

In this work we consider the ferromagnetic phase, in which there are q degenerate ground states (with 
all sites having the same colour), and the elementary excitations are kinks corresponding to domain walls. 
The action of the scaling field theory of the system can be written as a perturbation of the conformal field 
theory (CFT) corresponding to the critical point: 
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oil) 
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where the couplings r and h are related to the lattice couplings 


r ~ 
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r\j 


T -T, 
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(2.3) 

(2.4) 


and X = (xi,X 2 ) are Euclidean time and space coordinates. The corresponding CFT |25) is defined for 
q < 4 and has the central charge |26) 

c(g) = 1 - , ^ , (2.5) 

where the parameter t is related to q via the relation 

_ TT (t — 1) , , 

and the thermal and magnetic helds e and a are identihed with spinless relevant primaries ^> 2 ,i and 
^(t-i)/ 2 ,{t+i )/2 of tho CFT and have scaling dimensions |26l 127] 


2 h('?) = = 


(2.7) 

2 /1^") = = 

(t-l)(t + 3) 

(2.8) 

8 t [t + 1) 


The thermal operator preserves the Sq symmetry while adding the magnetic held breaks the permutation 
symmetry according to Sq —>• Sq-i. 

In absence of the magnetic held the theory is integrable and the complete scattering theory is known 
|281l2^ fT]. In the ferromagnetic phase there exist q degenerate ground states and the elementary excitations 
are kinks interpolating between them, with their mass m related to the coupling r via the mass gap relation 

m 


where 


(2) _ _L 
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[r(i)r(i)]^/^ / r(-Dr(i) 
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= 0.1643033, 


(2.9) 

( 2 . 10 ) 


In a magnetic held h in the direction of the g-th colour the degeneracy between the q ground states is 
lifted. For h < 0 there is a single true vacuum and q — 1 metastable ones, while for h > 0 their roles are 
exchanged. The expectation values of the magnetic held operator in the direction of the y-th colour is in 
the a-th ground state of the zero held theory is written as 


vii) 

1*^7 (®)l 0«) = 

where can be calculated using the formulas in EH: 

= 1.3578383417---X (2.12) 

= 1.9382577836---X (2.13) 


As a result, the difference of the energy density between the false and the true vacua is given by |10| 

Ae = - (5eg ~ h (^(cTg)^ - a^q (2.14) 

The combination v^^'^q/ {q — 1) is denoted by and called the string tension; I3^^'>\h\ gives the slope of 
the linear potential induced by the magnetic held. 


3 










2.2 Mesonic and baryonic mass estimations 

Here we present the various known estimations for the masses in the confinement spectrum for the case 
of the Ising and the three state Potts model. 


2.2.1 Meson masses in the Ising model 

The first mass estimation is the one obtained by McCoy and Wu |2] 

= m(2 + (2.15) 


where —Zn is the nth zero of the Airy function and A is the dimensionless ratio 


w? 


(2.16) 


This solution can be derived from the quantum mechanical system of two kinks in a linear potential. The 
quantum mechanical system also allows solutions corresponding to the zeros for the derivative of the Airy 
function; however, the corresponding wave-functions are symmetric and so forbidden due to the fermionic 
nature of the kinks. 

The WKB mass spectrum can be obtained by solving the quantization condition 


sinh {2'dn) — 2'dn 


27r (n — 1/4) 
2m cosh (i9n) 


(2.17) 


It can be improved further by adding higher corrections in A: 


sinh (2-i?„) - 2-dn = 27r (n - 1/4) A-h ^ A^+^S'fc (??n) 
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The first term in this expansion is given in |U] and it is 
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The Bethe-Salpeter equation (with various improvements) leads to a low energy expansion of the form 



4m^ 


1 + 

k=l 


(2.19) 


with the parameter t = A^^^. Different approximations of the /ifc coefficients taking into account multi¬ 
quark corrections (such as quark mass renormalization and renormalization of the short range quark- 
antiquark interaction), and string tension renormalization can be found in |3l ITJ (H [8]. The low-energy 
expansion for the meson mass from the Bethe-Salpeter equation has the form 


4m^ 
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while the radiative corrections modify this expression according to 


m(f — 
m 


= + ^(4c2 - a2)t^ - + 0(t^°) 

6 4 


( 2 . 21 ) 
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where 

02 = 0.0710809... C2 =-0.003889 ... (2.22) 

are the leading order quark mass and string tension renormalization corrections computed in |6] and 


B2 = 0.8 


(2.23) 


is the leading order interaction renormalization correction obtained in UM- 

2.2.2 Mesonic states in the three state Potts model 

For this latter model, the Bethe-Salpeter has not been carried out (albeit the setting was established in 
|16)i. so we quote the linear potential quantum mechanics and the WKB results. 

h<0 

For this sign of the magnetic held there is a single stable vacuum and two metastable ones. The 
two-kink conhgurations are 


K3a{^l)K^3{^2) 

where 3 is the stable vacuum, while a = 1, 2 are metastable vacua. Due to the presence of this degree of 
freedom allowed by a both antisymmetric and symmetric solutions are allowed. From the simple quantum 
mechanical picture of two-kink conhgurations it is clear that in the sector of zero total momentum the 
charge conjugation (C) parity of kink-antikink and therefore also meson states is equal to their parity 
under spatial rehections. 

The spectrum predicted by quantum mechanics in the linear potential is 


= 'm{2 + + O 

= m{2 + + O (a*/’) 


where A is the dimensionless ratio 

_ /?(3)|h| 

2 

and —Zn is the n-th zero of the Airy function, while —z!^ is the n-th zero of its derivative. 
For this case, the WKB quantization is given by |16) 


(2.24) 

(2.25) 


(2.26) 
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27r 2 arctan -|- iA{2'dn) + 0{X) even 

/sinh(z7r/3-|-i?)^ 

\sinh(z7r/3 — '&)) 
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(2.27) 

(2.28) 

(2.29) 

(2.30) 


This includes effects of nontrivial kink-antikink scattering, and therefore despite being semiclassical it goes 
beyond the simple quantum mechanical result above. 
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h>0 

For this sign of the magnetic field there are two stable vacua and one metastable. The allowed two-kink 
configurations are 

KU'»i)Kzp{^2) 


Charged meson states correspond to a 7 ^ /3, and their spectrum is given by the WKB quantization 
condition |16| is given by 


sinh {2'dn) — 2 ??n 
A 





(2.31) 


Note that charged meson single-particle states do not satisfy periodic boundary conditions on a circle due 
to a 7 ^ /3, therefore they cannot be observed in TCSA. 

On the other hand, neutral meson states do not exist, as they easily decompose under the process 


where for a = 1,2 one has /3 = 2,1, respectively. This process is allowed by the Chim-Zamolodchikov 
kink scattering amplitudes |Tj; also note that the kinks mediating between the stable vacua 1,2 are not 
confined. 


2.2.3 Baryon masses in the three state Potts model 

For the case /i < 0 all the kinks are confined. As a result, one may have three-kink bound states of the 
form 


A'31 (i?i) Ari2 ('^2) Ar 23 (193) 

corresponding to baryons and antibaryons. Both of these particles have the same spectrum due to charge- 
conjugation symmetry, and can be modeled in the form of a quantum mechanical three-body system. The 
low energy estimates for the baryon masses were recently obtained by Rutkevich m with the result 

M± = (2.32) 

where the ± correspond to parity under space reflection, with the following numerical values of e for the 
first three states: 

= 4.602 = 5.912 e(]“ = 7.098 

ej- = 6.650 = 7.734 Cg = 8.753 ^ ^ 

There are no baryons for /i > 0 as the kinks between the two stable vacua are not confined. 

3 TCSA, RG and extrapolation 

3.1 TCSA for the Ising and Potts field theories 

3.1.1 Scaling Ising model 

The Hilbert space of any conformal field theory can be decomposed into products of irreducible represen¬ 
tations of the left and right moving Virasoro algebras, which can be specified by giving their left and right 
conformal weights as 

= (3.1) 
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and every such sector corresponds to a primary field ^ . For the Ising model with central charge c 
the full Hilbert space is 


0 ® 

9 ’ 


1 1 
2 ’ 2 


©5j_ j_ 

16 ’ 16 

where the sectors on the hrst line are even, the ones on the second line are odd. The Hamiltonian is 


H = + T dxe + h dxa 


= 1/2 

(3.2) 

(3.3) 


(3.4) 


where 

e = 4>i 1 a = <l>j_ j_ (3-5) 

2’2 16’16 

Note that for this model the values h and —h are physically equivalent since they are related by the Z 2 
symmetry of the conformal held theory. 

For the Ising model we used the following level cut-offs with the dimensions of the truncated Hilbert 
space indicated below: 


n 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

dim 

77 

127 

213 

338 

551 

840 

1330 

1994 

3023 

4476 


3.1.2 Scaling 3-state Potts model 

The scaling limit of Potts model at the critical point is a minimal conformal held theory with central 
charge 

4 

(3.6) 
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j. The Hilbert space of the Potts model is the D 4 modular invariant 

Td = Tdo © Td+ © Td— © Tdi 


where 
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The D/^ conformal held theory is invariant under the permutation group S 3 generated by two elements Z 
and C with the relations 

= 1 = 1 czc = z-^ 


which have the signatures 


sign Z = -|-1 sign C = —1 


(3.9) 


(3.10) 


The sectors in Tdo of (3.7) are invariant under the action of the permutation group S 3 , the ones in Td± 
form the two-dimensional irreducible representation, which is characterized by the following action of the 
generators: 


C\±) = ±|=f) 

Z\±) = cos |±) ©sin |=f) (3.11) 
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while those in T-Li transform according to the signature representation of S 3 . 
The Hamiltonian is given by 


H = + T J dxe + h j dxa (3.12) 

where 

£ = $22 a = J, (3.13) 

5’5 15’15 

In our considerations an important role will be played by charge conjugation parity: the even sector under 
C consists of Ho and 'H+, while the odd sector consists of Hi and H-, with the dimensions as a function 
of the level cut-off n given below: 


n 

6 

7 

8 

9 

10 

11 

dim, even sector 

634 

1210 

2426 

4437 

8258 

14545 

dim, odd sector 

816 

1572 

3039 

5592 

10121 

(17904) 


In the case of the odd sector the left/right descendent levels must be different for some of the helds 
the right in order to get spinless helds. In our convention, the truncation level is chosen to agree with 
the smaller of the descendent levels. In the extrapolations we used level cut-offs from 6 to 11 in the even 
sector and from 6 to 10 in the odd sector. 

3.2 TCSA numerics: conventions 

The conformal Hamiltonian on a hnite circle of circumference R can be written as 

4’h = |(ro + £o-(5) 

In our TCSA calculations we consider the zero momentum sectors for which Lq — Lq 
level cutoff on the spectrum of states according to 

Lo<n (3.15) 

with n a positive integer. 

The computations are all performed in units of the kink mass m of the h = 0 model, which means the 
volume is measured in dimensionless units r = mR and the dimensionless energy levels are e = E/m. The 
hnite volume energy levels are given as functions ej(r), with i indexing the different levels; by convention 
the vacuum is taken to correspond to i = 0. Since we are interested in the mass spectrum, the relevant 
quantities are the relative energy levels 

ej(r) = ei(r) - eo(r) (3.16) 

One can also introduce the scaling function dehned by 

Mr) = ^Mr) (3.17) 

In the conformal held theory, the scaling function is a constant and is given by the eigenvalue of the 
operator 

Lq + Lq — (3.18) 

corresponding to the given level. In the off-critical model the leading term of ej(r) for large r is given 
by the bulk energy and is therefore generally linear (the Ising case is an exception with a logarithmic 
contribution), and so the scaling functions dj(r) grow as r^. On the other hand, the bulk contribution 
is universal for all levels, and in the absence of certain ultraviolet divergences (which is the case for the 
models considered here, as none of the perturbing operators has conformal weight /i > 3/4 |2n)l. the 
relative energy levels ei{r) go to a constant for large r. 


(3.14) 

= 0 , and impose a 












3.3 Renormalization group improvement 

To leading order in the level cutoff n, the cutoff dependence of the TCSA can be canceled by allowing 
the couplings to run according to renormalization group equations derived from second order perturbation 
theory. For a Hamiltonian of the form 

R = ^ (Ao + To - — j + ^ Aa / dx<^a(x) (3.19) 

the leading RG equations are |2ni [22] 

Ac(n) - Ac(n - 1) = ^ Aa(n)Ab(n)C'^b^—(1 + 0(l/n)) (3.20) 

a,b 

~ _ 

(27r)^“2^“ 

dabc da hfy ha 


where are the CFT operator product expansion coefficients: 


^>a(2,2;)^6(0,0) 


^ab^cjO, 0 ) 


(3.21) 


In the perturbation theory calculation, however, large denominators may appear due to the dependence 
of the scaling functions. This can be compensated for by taking into account the universal part of all 
scaling functions, by modifying the RG equations following the prescription in |221 (23] 


Xc{n) - Ac(n - 1) = 


1 




2 n — do{r) 


^ Aa(n)Ab(n)Cgfc ^^^ ^ (1 + 0(l/n)) 


a,6 


(3.22) 


where the vacuum scaling function do{r) can be estimated by its TGSA value at the starting cutoff for the 
RG run. In a unitary field theory it can be argued that the vacuum scaling function is always negative, 
so this does not introduce any new singularities. In essence this modification means regrouping some 
potentially large 1/n corrections into the RG flow. 

Note that this prescription also gives a running coupling for the identity, which leads to an additive 
renormalization constant for all energy levels. Since we consider only relative energy levels, this contribu¬ 
tion can be omitted. 


3.4 Extrapolation 

The higher 1/n terms give state dependent corrections corresponding to non-local counter terms |2n | l22 | l^. 
Their construction is quite involved and they are not known in a fully analytic form yet. However, there 
is an efficient shortcut that is sufficient for the purposes of the present work. The leading terms already 
incorporated in the RG equations yield cut-off corrections after integration with the scaling form 

(3.23) 


corresponding to the occurrence of in the OPE 

In the Ising model, the leading exponents can be summarized as: 
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The exponent 0 corresponds to a logarithmic divergence in the ground state energy, which cancels 
from the relative energy levels. The exponent —7/4 also corresponds to ground state renormalization and 
therefore also cancels. The two other exponents are taken care of by the running couplings. Therefore 
it is only necessary to take into account the highest subleading 1/n corrections, which lead to a residual 
cut-off dependence of the form 



n 


where both subleading terms come from 1/n corrections to the eel term. 
In the Potts model, the leading exponents are 


(3.24) 


^a^b\^c 
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2 i) 
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ee 
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(where some fields contained in the OPE which give even higher exponents have been omitted) and the 
residual cut-off dependence is 




(3.25) 


Our prescription for the RG-TCSA is as follows. In units of m, r is just a fixed dimensionless number 
given by Therefore the TCSA has two dimensionless parameters, one of which is given by value of h 
in units of m, i.e. the ratio 

h = (3.26) 

and the dimensionless volume parameter r = mR. For any value of r and m the physical values of the 
perturbed OFT couplings are 


A. = - 


A. = 




(2vr: 

h 


A-2h. 


(3.27) 


^2—2ha 


m 


2-2ha 


(2vr) 


l-2h„ 


Taking these as initial conditions at re = oo, the couplings can be run according to the RG equations (3.22) 


to determine their value at the given cut-off re. In practice one can approximate the difference equations 
by substituting 

\c{n) - Ac(re - 1) ^ (3.28) 

and solving the resulting differential equations numerically. We remark that in all our calculations the 
couplings ran very little so this has practically no effect, but as a matter of principle it must be done before 
we proceed to extrapolation. Once the RG eliminated all the leading cut-off dependencies, the renormalized 
TGSA Hamiltonian can be numerically diagonalized and then the residual cut-off dependence eliminated 
by fitting (3.24) for the Ising and (3.25) for the Potts case. In the case of the Ising model it turns out 


that the residual cut-off dependence alternates in sign between odd and even cut-offs, so the data for even 
and odd values of re were fitted separately, as demonstrated in figure [TT| For the case of the Potts model 
no such alternation was observed, and the data could be reliably extrapolated including both even and 
odd values of the level cut-off re as illustrated in 13.21 We also remark that in the Potts case we took into 


account the two exponents indicated in (3.25); the 11/5 and 12/5 exponents are too close together, and 
their effect is to small compared to the leading 7/5 to include both in the fit. 
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(a) First excited state (b) Second excited state 



Figure 3.1: Extrapolation fits of the relative energy levels in the Ising model in the ferromagnetic phase 
with magnetic field h = 0.008 for the first four excited states at dimensionless volume mR = 10. 
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(a) First excited state (b) Second excited state 




(c) Third excited state 


(d) Forth excited state 


Figure 3.2: Extrapolation fits of the relative energy levels in the Potts model in the ferromagnetic phase 
with magnetic field h = —0.05 for the first four excited states at dimensionless volume mR = 10. 


12 






4 Results 


4.1 Testing ground: Ising model with magnetic field 


4.1.1 False vacuum 

From the TCSA data it is possible to evaluate the energy density of the false vacuum relative to the 
stable one. For smaller volume, the TCSA converges fast, but for greater volumes the efficiency of the 


extrapolation procedure is apparent. The theoretical predictions can be calculated using (2.14) where the 


renormalized string tension given in |B] has also been taken into account. However, at the present precision 


the two predictions cannot be distinguished. Our results are illustrated in Figure 4.1 


Note that the false vacuum level is not a continuous level in the volume, therefore its linear rise does 
not in fact contradict the statement that all relative energy levels ei{r) go to a constant for r —)• oo. The 
metastable states are seen as level avoidances in finite volume IS11I2]; this is best demonstrated by the 


false meson level in the Potts model, shown in subsection 4.2.2 For the false vacuum, however, the level 
avoidances are confined to small enough regions, so that the level can be glued together from the pieces. 
In addition, for small coupling the first level avoidance is at a much larger volume than shown in the plots. 


4.1.2 Meson masses 

The meson masses can be estimated by first extrapolating the levels separately for even and odd cut-offs, 
and taking the average of the two results. For a more precise measurement the numerical procedure was 


slightly modified by keeping only the n = 10 ... 15 data and using only the 1/n term in (3.24). The reason 
for this is that the meson level data do not allow fitting the 1 /v? term with a sufficient precision. 

Selecting a given extrapolated one-particle level, one then finds the volume where the level is the most 
flat. The data still contain exponential finite size effects, which can be suppressed by fitting the part of 
the extrapolated nth meson level just before the flat portion by an exponential function 


= mn + Ar, 


■,-Bnr 


(4.1) 


and taking nin as the estimated mass of the nth meson. The results are shown in Table 4.1 (we note that 


the data did not permit the volume extrapolation for h = 0.4, so the TCSA numbers quoted there are just 
the value of the level at the point where it is most flat). The theoretical prediction “Airy” is given by eqn. 
(2.15), while BS and WKB given by (2.17) and ( 2.21[ ), respectively; the last line (iWKB) corresponds to 
taking into account the first correction in (2.18). One can see that the prediction from simple quantum 
mechanics in a linear potential (“Airy”) is generally only precise to a percent level or even worse for large 
magnetic fields, while WKB is an order of magnitude better. The theoretically rather involved Bethe- 
Salpeter approach only improves on WKB for very low mass mesons, while for higher masses the WKB 
is generally better. It is also clear that for h > 0.2 the only working theoretical framework is the WKB 
method. 

The upshot is that for all practical purposes WKB can be taken as the most reliable description of the 
spectrum over all the parameter range: it gives an approximation within 10“^ relative precision. This is 
an important lesson given that there are no Bethe-Salpeter predictions available in the three-state Potts 
case yet; however, we can take the WKB as an accurate prediction for comparison with the three-state 
Potts meson data. 


4.2 Three-state Potts model 

For the three-state Potts model we only consider the domain h < 0, since according to the discussion 
in subsection |2.2.2 there are no meson one-particle levels for h > 0, and reading off meson masses from 
two-particle states is difficult both for numerical reasons (the spectrum is dense, so level identification 
is difficult) and for theoretical reasons (extraction of masses with any precision requires modeling the 
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(a) h = 0.008 



(b) h = 0.03 



(c) h = 0.05 

Figure 4.1: Energy of the false vacuum for different values of the magnetic field in the Ising model. Lines 
are the theoretical predictions: blue is for (2.14), while purple one takes into account the string tension 
renormalization (the two lines are almost indistinguishable in the graphs). Dotted lines for different cut-off 
data, red dots are extrapolated (average of even and odd cut-off extrapolations) 
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h 


mi 

m 2 

m 3 

0.0175 

TCSA 

2.303 

2.526 

2.71 

Airy 

2.3068 

2.5364 

2.7243 

BS 

2.3021 

2.5228 

2.7005 

WKB 

2.3000 

2.5223 

2.7003 

iWKB 

2.3011 

2.5225 

2.7004 

0.0250 

TCSA 

2.382 

2.659 

2.89 

Airy 

2.3891 

2.6803 

2.9187 

BS 

2.3816 

2.6590 

2.8813 

WKB 

2.3791 

2.6583 

2.8812 

iWKB 

2.3803 

2.6586 

2.8812 

0.0375 

TCSA 

2.497 

2.855 

3.15 

Airy 

2.5099 

2.8915 

3.2039 

BS 

2.4969 

2.8556 

3.1414 

WKB 

2.4941 

2.8552 

3.1418 

iWKB 

2.4954 

2.8553 

3.1417 

0.05 

TCSA 

2.597 

3.027 

3.37 

Airy 

2.6177 

3.0800 

3.4584 

BS 

2.5987 

3.0281 

3.3684 

WKB 

2.5958 

3.0283 

3.3701 

iWKB 

2.5971 

3.0282 

3.3698 

0.1 

TCSA 

2.933 

3.588 

4.11 

Airy 

2.9805 

3.7143 

4.3151 

BS 

2.9312 

3.5851 

4.0919 

WKB 

2.9317 

3.5946 

4.1130 

iWKB 

2.9320 

3.5933 

4.1114 

0.2 

TCSA 

3.447 

4.451 

5.22 

Airy 

3.5565 

4.7213 

5.6750 

BS 

3.4115 

4.3556 

5.0251 

WKB 

3.4474 

4.4508 

5.2257 

iWKB 

3.4426 

4.4455 

5.2209 

0.4 

TCSA 

4.22 

5.71 

6.85 

Airy 

4.4707 

6.3198 

7.8337 

BS 

3.8964 

4.8679 

4.9503 

WKB 

4.2293 

5.7277 

6.8706 

iWKB 

4.2074 

5.7102 

6.8559 


Table 4.1: Meson masses in the Ising model. Theoretical predictions are shown with 4 digits accuracy, 
while for TCSA we show the digits that can be reliably extracted (with the last digit having an estimated 
precision of order 1). 
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meson-meson scattering). In addition, as discnssed in snbsection 2.2.3 there are no baryon states for 
h > 0. 


4.2.1 False vacuum 


The relative energy of the false vacnnm against the volnme is shown in Fignre 4.2 
extrapolation is again very effective. 


Notice that the 


4.2.2 Meson masses 

In contrast to the Ising case here we can only nse the method of extracting the valne at the flattest portion 
of cnrve to estimate meson masses, dne to the presence of the “false meson” resonance corresponding to a 
kink-antikink bonnd state conhgnration starting and ending in one of the false vacna. The “wavy” featnre 


these resonance plateans iHiiia introdnce in the spectrnm prevent application of the exponential ht ( |4.1[ ) 
to eliminate finite size effects in the meson mass. This effect can be seen in the plot |4.3[ which also 
demonstrates the efficiency of the nnmerical extrapolation procednre. The mesonic spectrnm against the 
absolnte valnes of the magnetic held can be seen on hgnrejT^ The deviations between the WKB prediction 
and the nnmerically determined masses are typically of the order of a few times 10“^, except in a few cases 
when a larger deviation of order 10“^ is observed. These are cases when the hattest portion of the meson 
level contains a level crossing with the false vacnnm, which makes the trnncation level extrapolation less 
precise. 


4.2.3 Baryon masses 

The baryonic states are in the higher part of the spectrnm. As a resnlt, becanse of the many level crossings 
the state mnst be carefnlly identihed for each valne of the volnme and cnt-off in order to carry ont the 
extrapolation. The masses are extracted as the valne of the extrapolated energy levels at its flattest point. 
As noted in snbsection 2.2.3| baryons and antibaryons have the same spectra in inhnite volnme. In finite 
volnme the eigenstates are the charge conjngation (C) even and odd combinations, as can be seen from 
the resnlts shown in Fignre |4.5| The deviations between the theoretical prediction and the nnmerically 
determined masses are typically of the order of a few times 10“^, except in a few cases when a larger 
deviation of order a few times 10“^ is observed, which in this case is mostly dne to difHcnlties of locating 
the level in the dense part of the spectrnm. 


5 Conclusions 

In this work we investigated confinement in the g-state scaling Potts field theory, for the cases q = 2 
(Ising) and q = 3 (three-state Potts). While these phenomena in the Ising model have been investigated 
in nnmerons works IS 0 m El [9] since the seminal paper by McCoy and Wn j2], resnlting in a very de¬ 
tailed nnderstanding of the meson spectrnm, in the case of the three-state Potts the theoretical predictions 
are more recent, especially regarding the baryon spectrnm. Onr method of choice was the renormaliza¬ 
tion gronp improved trnncated conformal space approach (RG-TCSA), a Hamiltonian trnncation method 
applied to pertnrbed conformal held theories, since there is no alternative for the three-state Potts model. 

Ising model was nsed both as a benchmark of the method, and as a way of comparing the effectiveness 
of theoretical predictions. Onr conclnsion was that semiclassical (WKB) qnantization was efficient over 
all the range of weak magnetic held, and therefore conld be taken as reference for the analysis of the Potts 
meson spectrnm. Indeed, we conld demonstrate very good agreement between the WKB and the meson 
spectra predicted in |16j . In addition, we compared the recent predictions the nnmerical resnlts to recent 
predictions for the baryon spectrnm and fonnd complete agreement. 

The present resnlts lead to the concinsions that from a qnantitative phenomenological viewpoint, the 
meson spectrnm of the the g-state scaling Potts held theory in weak magnetic held is described by WKB 
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(a) h = -0.005 



(b) h = -0.01 



(c) h = -0.02 


Figure 4.2: Energy of the false vacuum for different values of the magnetic held in the three state Potts 
model. Continuous lines are the theoretical predictions from (2.14), dashed lines are data with coupling 
constant renormalization before extrapolation for some values of the cut-off. Red dots are extrapolated 
data in the C-even sector, while black squares are extrapolated data from the C-odd sector. 
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Figure 4.3: Effect of the extrapolation for the hrst three mesonic state in the C-even and odd sectors with 
magnetic held h = —0.02. Dashed lines for data with running coupling for some values of the cut off, 
while large dots with dotted lines are the extrapolated data. We also plotted the energy of the even sector 
false vacuum (both before and after extrapolation, the latter marked with green) shifted up by 2 m to 
demonstrate that the wavy feature corresponds to the “false meson” resonance (meson conhguration over 
the false vacuum). 
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1^1 

Figure 4.4: Meson masses against the magnetic field in the three state Potts model. Black dots show 
TCSA results from the C-even sector while red is for C-odd. The green and blue lines show the WKB 
predicted for parity even and odd states. As expected, for mesons the two parities coincide. 
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Figure 4.5: The three lowest baryon masses against the magnetic held in the three state Potts model. 
Black dots show TCSA results from the C-even sector while red is for C-odd. The green and blue lines 
show the predictions ( 2.32[2.33 ) for space parity even and odd states, respectively. It can be seen that 
for each spatial parity there are two states, corresponding to the C-even/odd combinations of the baryon 
with its antiparticle. 
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to a very high precision, while the baryon spectrum can be efficiently modeled by the three-body quantum 
mechanical model introduced in m These findings are expected to be relevant in future investigations 
of the Potts field theory as a description of statistical systems. 
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